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The glass transition of mesoscopic charged particles in two-dimensional confinement is studied 
by mode-coupling theory. We consider two types of effective interactions between the particles, 
corresponding to two different models for the distribution of surrounding ions that are integrated 
out in coarse-grained descriptions. In the first model, a planar monolayer of charged particles is 
immersed in an unbounded isotropic bath of ions, giving rise to an isotropically screened Debye- 
Hiickel- (Yukawa-) type effective interaction. The second, experimentally more relevant system 
is a monolayer of negatively charged particles that levitate atop a fiat horizontal electrode, as 
frequently encountered in laboratory experiments with complex (dusty) plasmas. A steady plasma 
current towards the electrode gives rise to an anisotropic effective interaction potential between the 
particles, with an algebraically long-ranged in-plane decay. In a comprehensive parameter scan that 
covers the typical range of experimentally accessible plasma conditions, we calculate and compare 
the mode-coupling predictions for the glass transition in both kinds of systems. 

PACS numbers: 64.70.Q-, 66.30.jj, 64.70.ph, 64.70.pe 


I. INTRODUCTION 

Two-dimensional (2D) configurations of mesoscopic 
charged particles can be observed in various kinds of ex¬ 
periments [I| , including colloidal suspensions confined to 
interfaces or between plates ii, or negatively charged 
dust particles levitating in the weakly ionized plasma 
sheath atop and parallel to a flat horizontal electrode 
[^. In coarse-grained descriptions one is interested in 
the charged particle’s dynamics and phase behavior with¬ 
out taking explicit account of the surrounding electrons 
and ions that ensure overall charge-neutrality of the sys¬ 
tem. In this article, we employ mode-coupling theory 
(MCT) to study vitrification in two kinds of confined, 
monodisperse charged-particle model systems. The first 
is the traditional two-parametric model of confined par¬ 
ticles that interact via screened Coulomb (Yukawa) pair- 
potentials, and the second is a more realistic, three- 
parametric model for a monolayer of negatively charged 
particles embedded in a flowing plasma. 

The simple Yukawa model has been widely used in the 
description of dusty plasmas (see Refs. [l|, [a 0). It is 
capable of describing the effective pair-potential between 
charged particles rather accurately around the most com¬ 
mon (mean geometric) nearest neighbor distance 0, Q . 
Nevertheless, the Yukawa model is not justified in many 
of the common laboratory experiments with 2D confine¬ 
ment, due to a highly anisotropic distribution of ions. 
In the common case of dusty plasmas, levitating in a 
collisional plasma sheath atop an electrode in a radio 
frequency chamber Q, account has to be taken of the 
plasma current of ions towards the electrode and the 
corresponding anisotropic effective dust interaction po¬ 


tentials. A kinetic theory of the ion distributions and ef¬ 
fective dust grain interactions is appropriate in this case, 
and has been studied by different groups of researchers, 
under different assumptions on the plasma parameters 
The theory is based on the solution of the ki¬ 
netic equation for ions moving in the electrostatic field 
of the sheath. Different approximations used for the ion 
collision operator (describing the interaction with neu¬ 
tral gas) merely reflect different experimental regimes (in 
terms of the radio frequency discharge power and pres¬ 
sure) when the particular model is applicable. 

Among these kinetic models, the one published by 
Kompaneets et al. @ is based on a reasonable as¬ 
sumption of a mobility-limited ion drift in the sheath 
field (as opposed to rather unrealistic inertia-limited mo¬ 
tion) and employs a velocity-independent ion-neutral col¬ 
lisional cross-section which is logarithmica^ accurate for 
the dominant charge-exchange collisions The result¬ 
ing three-parametric potential is anisotropic in three di¬ 
mensions (3D); for charged particles confined to 2D, it 
exhibits an algebraically long-ranged r~^ decay. This 
model is expected to provide a realistic description of in¬ 
teractions in ground-based dusty plasma laboratory ex¬ 
periments El- 

Our results, reported in the present paper, predict 
qualitatively similar liquid-glass transition curves for 
monolayers with Yukawa-like and Kompaneets-like pair 
potentials. However, we find that a glass transition in a 
dusty plasma monolayer may be qualitatively misinter¬ 
preted if Yukawa-like interactions are assumed: An ap¬ 
parently re-entrant liquid-glass-liquid state sequence is 
found in the parameter space of the Yukawa potentials 
that at distances close to the mean geometric distance 
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best fit the potential derived from the kinetic theory. 
This apparent re-entrant state sequence is merely an arti¬ 
fact that arises when one attempts to describe the system 
in terms of the inappropriate Yukawa potential param¬ 
eters, and it disappears when the more realistic kinetic 
potentials are assumed, and the corresponding dimen¬ 
sionless parameters are used in plotting the transition 
diagram. 

The article is organized as follows: In Sec. |ll] we dis¬ 
cuss the two model systems of charged particle monolay¬ 
ers with Yukawa and Kompaneets pair potentials. Sec¬ 
tion mil provides a brief summary of the MCT equa¬ 
tions and their only input, the 2D static structure fac¬ 
tors, which are computed in the approximate T/2-HNC 
scheme. Our results are presented in Sec. El preceding 
our finalizing conclusions in Sec. [V] 


II. THE TWO MODEL SYSTEMS 

Both model systems that are described in the follow¬ 
ing two subsections contain mesoscopic charged particles 
confined to a 2D plane. The charged particles’ diam¬ 
eter is in the order of microns. Surrounding ions are 
only implicitly accounted for, through their influence on 
the effective pair-potential between the confined, charged 
particles. In the thermodynamic limit, both the number, 
N, of particles and the area, of the confining plane 
diverge to infinity at a fixed value of the areal particle 
number density n = N/L'^. 


A. Yukawa monolayer 




FIG. 1. Edge-on schematic of a Yukawa monolayer. Charged 
particles (filled circles) are confined to a plane, while op¬ 
positely charged ions are free to move in the surrounding, 
unbounded 3D space. The mean ion density is color-coded. 
Typical in-plane nearest neighbor distances are similar to the 
mean geometric distance and of the same order of mag¬ 

nitude as the Yukawa screening length Ay. Particle sepa¬ 
rations greatly exceed the particle diameter. The effective 
particle interactions are quantihed by the two dimensionless 
parameters Py and Ky = l/{XY\/n). 

The Yukawa monolayer model implicitly assumes ther¬ 
modynamic equilibrium statistics of ions, as schemati¬ 
cally depicted in Fig. [TJ Unlike the two-dimensionally 
confined, mesoscopic charged particles, ions are free to 


move in 3D space in the absence of external forces. Under 
these conditions the effective interaction potential Uy (x) 
between charged particles at sufficiently large mutual dis¬ 
tance follows the screened Coulomb (Yukawa)-type form 


Uy{x) exp(-Kya;) 

= ^ 

where ks is Boltzmann’s constant, T is the absolute tem¬ 
perature, and X = ry/n is the particle center-to-center 
distance in units of the mean geometric distance . 

The Yukawa potential in Eq. o is characterized by 
the two dimensionless parameters Ty and Ky'. The cou¬ 
pling parameter Ty = Q\^Jnl(An^tk^T') quantifies the 
interaction strength in terms of the charged particle’s ef¬ 
fective Yukawa charge Qy (which is tYpically less than 
the bare electric charge of the particles [ 13 , l^), and the 
dielectric permittivity e of the embedding medium. In 
case of dusty plasmas, e is equal to the dielectric per¬ 
mittivity of vacuum, cq, for all purposes of the present 
article in which we adhere to SI units. The screening 
parameter, Ky = l/(Ayy/n), is the normalized inverse 
of the Debye screening length Ay, which depends on 
the ion population. In an embedding plasma that con¬ 
sists of neutral partic les and univalent positive ions only, 
Ay = y/eksT/{e^rii) is the Debye length in terms of the 
proton elementary charge e, and the unperturbed (3D) 
ion number density rii of the ions far from the charged 
particle’s confining plane. 

The Yukawa model in two dimensions is best real¬ 
ized experimentally for charged colloids which are con¬ 
fined between two highly charged glass plates imm. 
There, the screening is caused by the microions between 
the plates 0 and it can be tuned by adding salt. The ex¬ 
perimentally observed freezing phase sequence has been 
found to agree with the theoretical predictions assuming 
a 2D Yukawa interaction 12211 . 


B. Kompaneets monolayer 

The second class of systems studied in this article is 
schematically depicted in Fig. [2] A radio frequency dis¬ 
charge chamber contains a weakly ionized plasma (of neu¬ 
tral gas particles, electrons, and ions), and negatively 
charged dust particles are levitating atop an electrode 
on the bottom of the chamber. Confinement of the dust 
particles to a well-defined 2D layer is achieved by a force 
balance between gravitation and electrostatic repulsion. 
Unlike the particles in the spatially unbounded Yukawa 
system, the ions in the radio frequency chamber exhibit 
a highly non-equilibrium steady state with a non-zero 
plasma current towards the electrode, where positive ions 
are adsorbed. Attraction between dust particles and 
ions causes downstream focusing of ions in the so-called 
plasma wake region. As a consequence, every dust par¬ 
ticle trails a positive space-charge in the downstream di- 
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rection, which causes the effective pair-potential between 
charged dust particles to be anisotropic in 3D. 



Electrode 


FIG. 2. Edge-on schematic of a dusty plasma monolayer. 
Negatively charged dust particles (filled circles) levitate in a 
well-defined 2D layer above an electrode in a radio frequency 
discharge plasma chamber, at a height where gravity is bal¬ 
anced by the vertical electrostatic force. The mean distri¬ 
bution of ions is color-coded. Three characteristic ion tra¬ 
jectories are sketched by arrows signed with -I-. Subsequent 
collisions between ions and neutral particles are separated on 
average by the ion-neutral mean free path 1 . Ions are focused 
in the downstream direction below the dust particles, giv¬ 
ing rise to positive space-charges in the plasma wake region. 
The effective dust particle interactions are quantihed by the 
three dimensionless parameters Tk, = 1/{\K\/n), and 
C = Xk/1. 


Kompaneets et al. Q have presented a self-consistent 
steady-state solution for the effective pair potential be¬ 
tween the dust particles, taking into account the exter¬ 
nal electric field E towards the electrode, and the colli¬ 
sions between ions and electrically neutral particles in the 
plasma. The resulting effective particle pair-potential, 
obtained under the assumptions of the mobility-limited 
ion drift in the field E, velocity-independent ion-neutral 
scattering cross section, and further assumptions that are 
outlined in the original reference, has been derived and 
described comprehensively in Ref. Q. We will refer to 
this kinetic pair-potential as the Kompaneets pair po¬ 
tential. For particles that are perfectly confined to a 
plane perpendicular to the plasma current, the in-plane 
Kompaneets potential Uk(x) is given by 


Uk{x) 

kBT 


Tk 


2C^iy 

TT 



dt 

l + C-^Y{t) 


X Ko 



1 + C-^Yit) ) ’ 


( 2 ) 


where Kq is the zeroth-order modified Bessel function of 
the second kind , and the two auxiliary functions X {t) 


and Y (t) are defined as 

Xit) = 1 - VTTii, 

_ 2y/TTii da 1 (3) 

Tt X [l + ff(l-a2)]2 “ + 

In Eq. m, the prefactor Tk = Q\^Jnl quan¬ 
tifies the interaction strength in terms of the effective 
charge Qk-, and the screening pa rameter is defined as 
kk = where \k = \/eEl/{e'^ni) is a field- 

induced screening length. In addition, the Kompaneets 
potential depends on the collision parameter C = Xk/1, 
where I is the mean free path between two consecutive 
collisions of an ion and neutral gas particles (“ion-neutral 
mean free path”, for short). 

For close-contact configurations (r <C C*ifwhere 
1/3 < s < 1, depending on the magnitude of C 0), the 
Kompaneets potential tends to the bare Coulomb poten¬ 
tial: 


X < CV«if : 


Uk{x) Tk 
k^T X 


(4) 


Hence, the Coulomb potential is recovered at all distances 
X in the limit C —>■ oo, corresponding to a very large 
field E, or a very small ion mean free path I or/and ion 
density rii. For large particle separations and finite values 
of C, the Kompaneets potential reduces to its in-plane 
asymptotic form 


Uk{x) 


kBT 


+ (5) 


The leading order asymptotic form of the anisotropic out- 
of-plane electrostatic potential is proportional to a;“^, 
and is given in Eq. (8) of Ref. 0 (in Gaussian units). 

In typical dusty plasma experiments the effective inter¬ 
action potential can be measured for particle distances 
X ~ 1 (ie., close to the mean geometric distance) W 
particle video tracking 0. It has been shown in Ref. 0, 
that the pair-potential in the experimentally directly ac¬ 
cessible narrow range of particle separations can be fitted 
equally well by the Yukawa as well as the Kompaneets 
form. However, one should expect that the qualitative 
differences between the Yukawa and Kompaneets poten¬ 
tials, most particularly in their long-ranged asymptotic 
forms, can have a considerable influence on collective dy¬ 
namics da and phase transitions. 


III. MODE COUPLING THEORY 

The glassy state is characterized by liquid-like static 
pair correlations without long range order, and a 
non-zero value of the non-ergodicity parameter fq = 
limt^oo (/)q(t), which is the long time limit of the 
wavenumber- and time-dependent autocorrelation func¬ 
tion 4>q{t) of the number density. The parameter fq is 
also called the form factor or the Debye-Waller factor. In 
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contrast to the glassy state, the liquid state is character¬ 
ized by a vanishing non-ergodicity parameter, fg = 0, for 
all wavenumbers q. In MCT, fg is calculated as [^, 1^ 

= (6) 

where in 2D [2^ 






(fk SkSp (q • k Cfe -I- q 


Sa f ,0, ^ ^ , , \2/,r 

P Cp) Jkjp, 

(7) 

with p = q — k. 

The static structure factor Sg and direct correlation 
function Cg = 1 — 1/Sg are the only input to the MCT 
equations, conveying information about the particle in¬ 
teractions. Note that the number density n does not 
explicitly enter into Eq. 0, since all lengths and wave 
vectors are expressed in units of l/\/n and y/n, respec¬ 
tively: In our notation the wave vector q is the dimen¬ 
sionless Fourier conjugate variable to the dimensionless 
distance vector x — r^Jn. 

The Lamb-Mossbauer factor /® = limt_,.oo which 

is the long-time limit of the wavenumber- and time- 
dependent, Fourier transformed tagged particle position 
autocorrelation function ipgit), is calculated in MCT ac¬ 


cording to [271 


l-/« 


= J^q[f,n 


( 8 ) 


where 


Sk{q ■ kfclfkf; 


(9) 


and, once again, p = q — k. 

We evaluate the integrals in Eqs. 0 and 0 numer¬ 
ically and solve Eqs. 0 an^8]) iteratively with itera¬ 
tion seeds /q° = /|° = 1 [ls|. To evaluate the inte¬ 
grals numerically we use N = 200 equidistant grid points 
with spacing Aq = 0.2, minimal wavenumber qmin = 0.1 
and maximal wavenumber (/max = 39.9. Test calculations 
with N = 500 grid points allow us to estimate the numer¬ 
ical error due to integral discretization, which is around 
5% in the glass transition temperatures. The static struc¬ 
ture factor is obtained from the (Fourier transformed) 
solution of the T/2-HNC integral equation [2^ 


7(a;) = J ( 


X C X — X 


exp < — 


2U{x') 

ksT 


- 1 


( 10 ) 

for an isotropic 2D fluid in terms of the indirect and 
direct correlation functions 'yix) and c{x) [s^. Equa¬ 
tion m is solved by means of a numerical spectral solver 
for liquid integral equations in an arbitrary number of 
spatial dimensions, which has been comprehensively de¬ 
scribed in Ref. [3l| , and which is based on methods that 
have been originally introduced in Refs. The 



FIG. 3. Static structure factors for a Yukawa monolayer with 
Fy = 100 and ky = 2.0 (lower three data sets) and a Kom- 
paneets monolayer with Tk = 300, kk = 2.0, and = 0.25 
(upper three datasets). Crosses and circles: Monte Carlo sim¬ 
ulation results. Dashed curves: HNC integral equation solu¬ 
tion. Solid curves: Solution of the T/2-HNC in Eq. (IIOII . 
The Kompaneets monolayer structure factors are shifted by 3 
units along the vertical axis for clarity. 


spectral solver operates on logarithmically spaced grids 
of wavenumbers and radii, providing a high-resolution 
structure factor that is mapped to the above mentioned 
equidistant wavenumber grid by quadratic interpolation. 

Note that Eq. m is a simple modification of the 
well-known hypernetted chain (HNC) integral equation 
pm, which is recovered when the term 2U{x') in the 
integrand is replaced by U{x'). Thus, the solution of the 
T/2-HNC equation coincides exactly with the solution 
of the HNC integral equation for a system in which the 
temperature has been scaled down by a factor of 1/2. 
In Ref. P, the MCT glass transition was studied for 
two-dimensional binary mixtures of aligned point-dipoles 
with a long-ranged repulsive pair potential that is pro¬ 
portional to the inverse cube of the particle separation, 
r~^. It was empirically found in Ref. P| that the T/2- 
HNC scheme predicts the static structure factors of the 
strongly repulsive 2D binary dipole mixtures with a sig¬ 
nificantly higher accuracy than the HNC scheme. In or¬ 
der to test the accuracy of the T/2-HNC scheme for the 
Yukawa- and Kompaneets monolayer systems, we have 
simulated 2D equilibrium liquids with strong repulsive 
pair potentials of both types, and compared the static 
structure factor from the simulation to the HNC and the 
T/2-HNC scheme solutions. Our results, shown in Fig.|3l 
underpin the good accuracy of the T/2-HNC and its 
supremacy over the HNC scheme. We have obtained the 
datasets represented by crosses and circles in Fig. [3] from 
Metropolis Monte Carlo (MC) simulations in the NLT- 
ensemble of constant particle number N, constant system 
area L^, and constant temperature T. A square sim¬ 
ulation box with periodic boundary conditions in both 
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Cartesian directions was used in our simulations, and we 
haven chosen the parameters Fy = 100 and Ky = 2.0 
for the Yukawa monolayer of IV = 10.000 particles and 
Tk — 300, kk = 2.0, and C = 0.25 for the Kompa- 
neets monolayer of IV = 12.000 particles. Both simulated 
systems are strongly coupled equilibrium liquids not far 
from the crystal-liquid transition point. In our MC sim¬ 
ulations, the direct particle interactions are truncated at 
a dimensionless cutoff radius of Xc = Tcsfn = 5 in case of 
Yukawa interactions, and at x^, = 12.5 in case of Kompa- 
neets interactions. For pair separations x > Xc, the pair- 
potential is set equal to zero in the simulations. Varying 
its numerical value, we have checked that the cutoff ra¬ 
dius is large enough and does not have a significant effect 
on the measured quantity Sq. 

Note from Fig. [3] that despite its improved accuracy 
in comparison to the standard HNC scheme, the T/2- 
HNC scheme still tends to underestimate the principal 
peak height in Sq. In addition note that we apply the 
T/2-HNC scheme in the following sections to systems 
at the liquid-glass transition, that is, beyond the equi¬ 
librium fluid regime for which the accuracy of the in¬ 
tegral equation scheme can be tested by comparison to 
crystallization-free simulations. Moreover, the approxi¬ 
mate T/2-HNC scheme is combined in the following with 
the approximate MCT equations. The combined uncer¬ 
tainty of the resulting glass transition lines cannot be 
easily estimated and, thus, the numerical values of the 
glass transition temperatures must be taken with some 
caution. Nevertheless, the dominating qualitative fea¬ 
tures of Sq are contained in the T/2-HNC solution, and 
the features of the glass transition curves can be expected 
to be at least qualitatively correct. 

It is important to note also that the T/2-HNC scheme 
is empirically justified only in case of strong enough par¬ 
ticle interactions. In the limit of vanishing interactions, 
Fy—^OorFif—^0, the T/2-HNC scheme predicts twice 
the correct asymptote c(r) —^ exp{—t/(r)/fcBT} — 1 for 
the direct correlation function (i.e., twice the Mayer func¬ 
tion) . A related issue is the wrong long-distance decay - 
the T/2-HNC scheme yields twice the correct expression 
limr_>.oo c(r) = —U(r)/kBT. This wrong long-distance 
decay is observed for all values of the potential prefactor. 


IV. RESULTS 
A. Glass transition diagrams 

The glass transition curves for the Yukawa monolayer, 
and for three different Kompaneets-monolayers with dif¬ 
ferent values of the collision parameter are shown in 
the transition diagram in Fig. IT] In the transition dia¬ 
gram, the screening parameters Ky and kk vary along 
the horizontal axis, and the coupling parameters Fy and 
Tk vary along the vertical axis. The data points in Fig.0] 
represent the lowermost values of Ty x for which fq as¬ 
sumes a non-zero value at given values of ky^k- Note 



FIG. 4. Glass transition curves in the (k, r)-plane, for the 
Yukawa potential (black curve with triangles) and three dif¬ 
ferent Kompaneets potentials with parameters ^ = 0.25 (solid 
curves with diamonds), ^ = 0.375 (solid curve with squares) 
and ^ = 0.5 (solid curve with circles). The dashed curve is 
the 2D Yukawa freezing line from Ref. . 


that at the glass transition, fq > 0 for all finite values 
of q smaller than the q cutoff. In our implementation, 
we have tested at g = 3.9 to identify the glass transi¬ 
tion points. In the One-Component-Plasma (OCP) limit 
Ky^K 0, both the Yukawa and the Kompaneets po¬ 
tential reduce to the unscreened Coulomb potential [see 
Eqs. (dl and (|4])], and the glass transition curves close in 
on the T/2-HNC-MCT approximation for the OCP glass 
transition point, Fy = F/f = 138.5. 

While the glass transition curves are qualitatively sim¬ 
ilar for the Yukawa and the Kompaneets systems, the 
transition occurs at higher values of the coupling param¬ 
eter in case of the Kompaneets monolayer. For decreas¬ 
ing values of the parameter the differences between 
the Yukawa and Kompaneets glass transition curves are 
increasing. Such trend is not surprising: As we discuss in 
the next subsection (see also Fig. ED, the deviation of the 
Kompaneets potential from the Yukawa-like form dras¬ 
tically increases as f decreases. On the other hand, in 
the limit f ^ oo the Kompaneets potential tends to the 
Coulomb form, so the Kompaneets glass transition curve 
in this case would be a horizontal line in the transition 
diagram of height F^ = 138.5. 

For 3D Yukawa systems it has been found that the 
glass transition and crystallization (freezing) lines are 
approximately parallel in the (Ky,Fy)-plane [3^ . The 
same similarity between the glass transition and crys¬ 
tallization lines is found for the 2D Yukawa monolayer 
in Fig. m where we plot the 2D Yukawa freezing line 
reproduced from Ref. (dashed curve), and the T/2- 
HNC-MCT 2D Yukawa glas s transition line (black curve 
with triangles). In Ref. |38l |. the crystallization line was 
obtained from simulations, and it was approximated by 
the inverse polynomial F = F*/(l -|- / 2 K^ -|- /sk^ -I- / 4 k‘*), 
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with r* = 73.9 = 131/y?, /2 = -0.1235, /a = 0.0248 
and /4 = —0.0014. Note that the 2D ion-sphere radius 
0=1/(also called Wigner-Seitz radius) was used as 
a unit of length in Ref , instead of the mean geometric 
distance 1/v^ utilized in the present paper. Therefore, 
one has to take account of a l/-\/7r prefactor difference 
in the definitions of the Yukawa coupling parameter and 
the inverse Yukawa screening parameter. 


B. Potentials and structure factors at the glass 
transition 

In Fig. [S]we plot the Yukawa potential and three dif¬ 
ferent in-plane Kompaneets potentials for different val¬ 
ues of C, all at the glass transition for ny = kk = 2.0. 
The full set of parameters, including the glass transi¬ 
tion values of rY and Ti^, is provided in the figure cap¬ 
tion. Note that the potentials in Fig. [5] are multiplied 
by their argument x = ry/n, to expose the differences. 
The curves corresponding to Kompaneets pair-potentials 
(with r~^ asymptotics) therefore decay proportionally to 
x~^ for large values of x. The inset of Fig. [5] features the 
T/2-HNC static structure factors Sq, corresponding to 
the four different potentials plotted in the figure’s main 
panel. Despite the pronounced differences between the 
four potentials (in particular around the most frequently 
sampled mean geometric distance x = 1), all four func¬ 
tions Sq at the glass transition are indistinguishable on 
the scale of the figure inset. The principal peak heights 
of the four structure factors differ only slightly in their 
values. 


C. A fallacious re-entrant state seqnence 

As pointed out in Sec. im the 2D Yukawa model 
with its many simplifying assumptions is merely a toy 
model for experimentally observable monolayers of meso¬ 
scopic charged particles. For the important class of 
ground-based dusty plasma experiments the kinetic pair- 
potentials are far more realistic, and among them the 
Kompaneets pair potential stands out with its realistic 
model assumptions. In this section, we allude to the 
possible consequences of an over-simplified interpretation 
of charged particle monolayers in terms of the Yukawa 
model. We show that the liquid-glass transition of a 
system with Kompaneets-like pair potential appears as 
a non-monotonic curve (corresponding to liquid-glass- 
liquid state re-entrance) when it is plotted in terms of the 
inappropriate parameters of the Yukawa potentials that 
represent a best fit to the actual (Kompaneets) potential 
around the mean geometric distance x = r^Jn = 1. 

In Fig. [B] we plot the Yukawa glass transition curve that 
is also shown in Fig. |4] (black curves with triangles). The 
2D Yukawa freezing line, reproduced from Ref. [38|, is 
also shown (dashed line) to allow a better comparison to 
the glass transition line than on the scale of Fig. 01 The 



FIG. 5. Effective pair-potentials for ky — kk = 2.0, for 
values of Fy and Fk at the liquid-glass transition point. 
Solid curve: Yukawa potential for Fy = 195.4. Dashed, dot- 
dashed and dot-dot-dashed curves: Kompaneets potentials 
for C = 0.25, 0.375 and 0.5, and F^ = 539.4, 406.8 and 345.0, 
respectively. All potentials are multiplied by their argument, 
X = r^/n, to expose the differences. The inset features the 
corresponding static structure factors Sq in T/2-HNC approx¬ 
imation. All four functions Sq are overlapping on the scale of 
the inset. The principal peak heights of the structure factors 
are Sq = 6.33 for the Yukawa system, and Sq = 6.26, 6.23 
and 6.19 for the Kompaneets systems with = 0.5, 0.375 and 
0.25, respectively. 


curve with open squares in Fig.[B]is generated as follows: 
For given values of the two Yukawa parameters Ky and 
Fy, we calculate the Kompaneets potential that fits best 
to the Yukawa potential in the distance range 0.7 < a; < 3 
which is most frequently sampled by the particles Q. 
The fit is conducted as follows: For given values of I and 
n, which yields the combination Qkk = {l\/n)~^ (~ 0.354 
for the example shown in the figure), we tune the two re¬ 
maining, independent Kompaneets parameters kk and 
Fif; an optimal fit is achieved by minimizing the square 
deviation jQ^dx[Uy{x) — Uk^x)]^ between the two po¬ 
tentials. We then calculate Sq for the best-fitting Kom¬ 
paneets potential in the T/2-HNC scheme, and use it as 
the input to the MCT equations (O and ([7]) for /,. If 
fq = 0, the system is classified as liquid, and if fq > 0, it 
is classified to be in the glassy state. We repeat the full 
procedure for various Yukawa parameters Ky and Fy, 
which are tuned by interval bisection, until we find for 
each Ky the smallest (critical) value of Fy at which the 
best-fitting Kompaneets system vitrifies. 

Thus, the curve with open squares in Fig. |6] is the 
glass transition curve of a dusty plasma monolayer with 
Kompaneets-like interactions, as it would appear when 
plotted in terms of the dimensionless parameters Ky and 
Fy of the Yukawa potentials that best fit the actual Kom¬ 
paneets potential around the mean geometric distance, 
where the potential is directly accessible Q. Therefore, 
if one observes vitrification in a dusty plasma monolayer 
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and assumes Yukawa-like interactions in the experiment 
analysis, the transition behavior may be misinterpreted 
as a re-entrant liquid-glass-liquid state sequence, while 
the transition diagram in terms of the three relevant 
Kompaneets potential parameters does not exhibit any 
re-entrance (see Fig. 01). 



FIG. 6. Kompaneets and Ynkawa glass transitions in the 
Ynkawa parameter plane. The cnrve with triangles: Ynkawa 
glass transition curve in the Yukawa screening- and coupling- 
parameter (fty, ry)-plane. The curve with squares is the glass 
transition curve for Kompaneets pair-potentials which have 
been optimally fitted to the corresponding Yukawa potential 
in the region 0.7 < a; < 3, by pointwise tuning of the Kom¬ 
paneets screening parameter kk and coupling parameter T k ■ 
The parameters I = 2.3 mm and n = 1.5 mm“^ are held fixed 
for the Kompaneets potential. The dashed curve is the 2D 
Yukawa freezing line from Ref. [^ . 


D. Non-ergodicity parameters 

We turn our attention now to the g-dependent Debye- 
Waller and Lamb-Mossbauer factors at the glass transi¬ 
tion, which are plotted in Fig. [T] It is observed that fg 
approaches zero in the limit g —^ oo. In the opposite limit 
g —>■ 0, the functions fg for Yukawa- and Kompaneets 
potentials with a small value of the screening parameter 
ky = Kk = 0.05 assume very small but non-zero val¬ 
ues, and for finite wavenumbers g, all plotted functions 
fg deviate clearly from zero. Finite wavelength density 
modulations, cannot relax in the glassy state since this 
would require a collective rearrangement of particles on 
the length scale of some nearest neighbor cage diameters. 
Observing Fig.[7]and the inset of Fig.[5l one can see that 
the most resilient density modulations (corresponding to 
the principal maximum in fg) are for g Ri 27r, that is, at 
the wavenumber that corresponds to the static structure 
factor principal maximum, and to the Fourier conjugate 
of the nearest neighbor (mean geometric) distance. Very 
long wavelength (g —> 0) density modulations cannot re¬ 
lax in the glassy state in general, as indicated by the finite 



0 10 20 30 40 

q 


FIG. 7. The Lamb-Mossbauer factors fg (monotonically 
decaying as functions of q) and the form factors fg (non¬ 
monotonic functions of q) for various Yukawa monolayers 
(symbols) and Kompaneets monolayers (curves) at their re¬ 
spective MCT glass transition points, with potential param¬ 
eters as indicated in the legend. 


values of fg^o in Fig. 0 This is due to the finite isother¬ 
mal compressibility, Sg^o, of the system. Note here in 
particular that the Debye-Waller factor of hard-sphere 
system remains finite as g approaches zero [2^, which is 
in line with the rather large compressibility of such hard- 
sphere system. However, in the OOP limits ky ^ 0 and 
Kk —>■ 0, in which both the Yukawa and the Kompaneets 
potential reduce to the Coulomb potential, the isother¬ 
mal osmotic compressibility coefficient vanishes [10, SO] , 
which corresponds to an infinite thermodynamic driving 
force for the leveling of long-wavelength density modula¬ 
tions. Therefore, fg-yo = 0 in the OCP limit. 

The Yukawa- and Kompaneets system are indistin¬ 
guishable in the OCP limit. This facilitates computa¬ 
tion of the OCP-limiting behavior in terms of a small-q 
asymptotic expansion of the Yukawa monolayer Debye- 
Waller factor. As demonstrated in the Appendix, the 
T/2-HNC solution for the direct correlation function of 
a Yukawa monolayer can be approximated as 


Cg 


2UY{q) 

keT 


dTrFy 

\/ky 


for g -I- gt > Ky + 9^, 


( 11 ) 


that is, when both the wavenumber g and the screen¬ 
ing parameter ky are small and within a certain ratio 
of each other. In Eq. dm), qt is a dimensionless non¬ 
negative threshold wavenumber with a typical value of 
qt 0.1. The corresponding small-g, small-Ky form of J^g 
is obtained from a functional Taylor expansion Hi, 113, 
resulting in 


J^g = {a + Pq^ .. .)Sg 


( 12 ) 










































where 5, = (1 — Cq) 

“ = ^ y dk kSl (^4 + kckc'k + fk ( 13 ) 

and 

P=^ [ dk ksl (c'4 + ^k^cf + h44 

8.J 32 2 

+Ckcl + -kckc'l' + -^k'^4c'l'j fl- 

Considering only the leading order of the approximation, 
this translates into the small-g, small-Ky limiting behav¬ 
ior of the Yukawa monolayer Debye-Waller factor, 


analysis: While the glass transition line is a monotonic 
function in terms of the three relevant, dimensionless 
Kompaneets pair-potential parameters, it appears to be 
non-monotonic corresponding to a fallacious liquid-glass- 
liquid re-entrance when the pair interactions are misin¬ 
terpreted as Yukawa-type interactions. 

A promising task for future research would be the gen¬ 
eralization of our results to binary systems, in order to 
understand the different glass types in mixtures. Since 
the crystalline states in such systems are pretty complex 
(dljl , the glass transition scenarios are also expected to be 
much more complex than in the monodisperse system. 

ACKNOWLEDGEMENTS 


fq = a 


1 -I- a -I- 


AttVy 


for q + qt'> Ky + q^, 


(15) 

in T/2-HNC approximation. For finite Ky, the function 
fq in Eq. (dS]) assumes a positive value for q = 0, and 
increases oc q^ when q —>■ 0. Only in the OCP limit 
Ky = 0, the function fq in Eq. m vanishes for q = 0, 
and increases initially as oc q. In a broad scale the fq 
asymptotic for Ky = 0.05 and Fy = 138.7 is almost 
linear. 

Note here that the small-q limiting OCP Debye-Waller 
factor is qualitatively different in two and three dimen¬ 
sions. In 3D, the function fq vanishes in the OCP limit 
Ky = 0 as oc q^ [1^. In contrast to fq, the Lamb- 
Mossbauer factor /* in MCT approximation does not 
critically depend on the form of the pair potential, since 
Tq in Eq. ® and also the small-q limit of Tq do not 
depend on Sq. 


V. CONCLUSIONS 

We have calculated the liquid-glass transition bound¬ 
aries in the state diagram spanned by the screening pa¬ 
rameters and coupling parameters of 2D monolayers with 
Yukawa- and Kompaneets-like pair potentials, in T/2- 
HNC-MCT approximation. While both types of systems 
exhibit qualitatively similar glass transition curves, there 
is a quantitative difference in the vitrification tempera¬ 
ture, which decreases as a function of the collision pa¬ 
rameter f of the Kompaneets pair-potential. Both the 
Kompaneets- and Yukawa-monolayer reduce to a two- 
dimensionally confined OCP in the limit of infinite Ay 
and Xk- 

In contrast to the over-simplifying 2D Yukawa model, 
the kinetic pair-potentials, including in particular the 
Kompaneets pair-potential, provide far more accurate 
descriptions of the interactions between dust grains in 
typical ground-based complex plasma experiments. We 
have demonstrated that a glass transition in a dusty 
plasma monolayer is prone to a qualitative misinter¬ 
pretation if the simple Yukawa model is invoked in its 
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APPENDIX 



FIG. 8. Function |a;c*''’^(a:)| for Yukawa monolayers and vari¬ 
ous potential parameters, calculated in T/2-HNC approxima¬ 
tion, is bound from above by 2ry exp{—x/2}. 


Here we validate the small-q, small-Ky result for the 
T/2-HNC direct correlation function of a Yukawa mono- 
layer in Eq. CD- We begin by noting that the function 
c(x) can be split into the sum 




( 16 ) 


of a short-ranged part, and the asymptotic long- 

ranged part, —2ry exp(—Kya;)/a:. In Eq. (fTOl) . the pecu¬ 
liar long-ranged asymptotics c{x —>■ oo) = —2U{x)/kBT 
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of the T/2-HNC scheme solution has been taken into ac¬ 
count {c.f., our discussion at the end of Sec. Iml) . The 
direct correlation function in wavenumber-space is cal¬ 
culated as the isotropic 2D Fourier transform (Hankel 
transform) 


OO 

Cq= 27 t j dx xc{x)jQ{qx) (17) 

0 

OO OO 

= 27r y dx xc^‘^\x)jQ{qx) — J dx e~'^^^Jo{qx), 


Cl 


C2 


binations of Fy and ky, and even in the OCP limit 
ky = 0. This finding, combined with the upper bound 
I Jo(a;)| < min{l, ^y2/{^^x)} for the envelope of the Bessel 
function, allows us to compute an upper bound 

(18) 

for the modulus of the function ci{q). Solutions for all 
Hankel transforms occurring in such computation are 
listed in Ref. [i^. Noting that C 2 = dirFy/y^Ky -|- q"^, 
we conclude that Ci in Eq. (IT71) is negligible if the condi¬ 
tion 


|ci| < min < 127 rF 


y> 


SttF, 

7^ 


q + qt ^ Ky + q^ 


(19) 


where Jq denotes the Bessel function of the first kind 
with index 0. Note from Fig.[5]that the T/2-HNC scheme 
solution for the function |a;c^®^(a;)| is bound from above 
by the function 2rY exp{—a;/2} for all reasonable com- 


is fulfilled, where qt 0.1 is a threshold wavenumber. 
For all combinations of q and Ky that fulfill Eq. ra, 
the T/2-HNC solution for Cq is well approximated by 
Cq ~ -C2 = dTrFy/-I- 
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